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We have theoretically investigated a general function photonic crystals (GFPCs) with defect layer, 
and choose the line refractive index function for two mediums A and B, and analyze the effect of 
defect layer's position, refractive indexes and period numbers on the transmission intensity and the 
electric field distribution. We obtain some new characters that are different from the conventional 
PCs, which should be helpful in the design of photonic crystals. 
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1. Introduction 

Photonic crystals (PC) are a new kind of mater 
which facilitate the control of the light [1]. PC exliil 
Photonic Band Gaps (PBG) that forbids the radiat 
propagation in a specific range of frequencies [2-6] . 
PBG forbids the radiation propagation in a specific ra 
of frequencies. The existence of PBGs will lead to mi 
interesting phenomena, e.g., modification of spontane 
emission [7-9] and photon localization [10]. Thus nun 
ous applications of photonic crystals have been propo 
in improving the performance of optoelectronic and 
crowave devices such as high-efficiency semiconduc 
lasers, right emitting diodes, wave guides, optical filti 
high-Q resonators, antennas, frequency-selective surfj 
optical limiters and amplifiers [11-18]. In the past 
years has been developed an intensive effort to study < 
micro-fabricate PBG materials in one, two or three di- 
mensions [19-21]. In Refs. [22-25], we have proposed a 
general function photonic crystals (GFPCs), which re- 
fractive index is a arbitrary function of space position. 
Unlike conventional photonic crystals (PCs), which struc- 
ture grow from two materials, A and B, with different 
dielectric constants ea and e^, and have obtained some 
results different from the conventional photonic crystals. 
In the paper. We have studied the general function pho- 
tonic crystals (GFPCs) with defect layer, and choose the 
line refractive index function for two mediums A and B. 
We obtain some results: (1) When the position of de- 
fect layer move behind in the GFPCs, or the refractive 
indexes of defect layer increase, the transmission inten- 
sity maximum of defect model decreases. (2) When the 
period number of GFPCs increase, the transmission in- 
tensity of defect model increase, and the defect model's 
width of half height become narrow. (3) Both the de- 
fect layer and it's position have the effects on the electric 
field. (4) In the structure {BA)^ D{BA)^ , when the pe- 
riod number N increase, the relative intensity of electric 
field increase. (5) When the refractive indexes of defect 
layer {ud) increase, the relative intensity of electric field 
was enhanced. 
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FIG. 1: The motion path of light in the medium of refractive 
index n{z). 



2. The light motion equation in general func- 
tion photonic crystals 

For the general function photonic crystals, the medium 
refractive index is a periodic function of the space posi- 
tion, which can be written as n{z), n{x, z) and n(x, y, z) 
corresponding to one-dimensional, two-dimensional and 
three-dimensional function photonic crystals. In the fol- 
lowing, we shall deduce the light motion equations of the 
one-dimensional general function photonic crystals, i.e., 
the refractive index function is n = n{z), meanwhile mo- 
tion path is on xz plane. The incident light wave strikes 
plane interface point A, the curves AB and BC are the 
path of incident and reflected light respectively, and they 
are shown in FIG. 1. 

The light motion equation can be obtained by Fermat 
principle, it is 

6 [ n{z)ds = 0. (1) 

J A 

In the two-dimensional transmission space, the line ele- 
ment ds is 

ds = ^/{dxY -V [dzf = ^/^TPdx, (2) 
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FIG. 2: The light transmission and electric magnetic field 
distribution figure in FIG.l medium. 



where z = j^, then Eq. (1) becomes 

fB 

6 / n{z)^/TT{i)^dx = 

J A 

The Eq. (3) change into 



(3) 



dz 



dz 



At the two end points A and B, their variation is zero, 
i.e., Sz{A) = Sz{B) = 0. For arbitrary variation Sz, the 
Eq. (4) becomes 



dn{z) 
dz 
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simpliiy Eq. (5), we have 

dn{z) zdz 



n{z) 



(5) 



(6) 



The Eq. (6) is light motion equation in one-dimensional 
function photonic crystals. 

3. The transfer matrix of one-dimensional gen- 
eral function photonic crystals 

In this section, we should calculate the transfer matrix 
of one-dimensional general function photonic crystals. In 
fact, there is the reflection and refraction of light at a 
plane surface of two media with different dielectric prop- 
erties. The dynamic properties of the electric field and 
magnetic field are contained in the boundary conditions: 
normal components of D and B are continuous; tangen- 
tial components of E and H arc continuous. We consider 
the electric field perpendicular to the plane of incidence. 



and the coordinate system and symbols as shown in FIG. 
2. 

On the two sides of interface I, the tangential compo- 
nents of electric field E and magnetic field H are contin- 
uous, there are 



Eq = Ej = Eti + E'^2 

Ho=Hi = Hn cos 0^ - H'^^ cos 91. 



(7) 



On the two sides of interface II, the tangential compo- 
nents of electric field E and magnetic field H are contin- 
uous, and give 



Eji = E'j = Ei2 + Er2 

Hn = H'j= Hi2 cos 6»f - Hr2 cos 6»f , 



(8) 



the electric field En is 

En = Siioe*('=»"^+'=^^^U=o = i;tioe'^"(°)^^"^*^^, (9) 
and the electric field Ei2 is 

p 771 i(k'xB+k'z)i 

Where xa and Xb are x component coordinates corre- 
sponding to point A and point B. We should give the 

relation between Ei2 and En- By integrating the two 
sides of Eq. (6), we can obtain the coordinate compo- 
nent Xb of point B 



r^'^ dn{z) _ zdz 

Jn{0) n{z) 1 + i 



ko 



to get 



and 



dx 



dz 



(11) 



(12) 



(13) 



where fco = cot Of and kz — 
n{z) > n(0) ahiOl. and the coordinate Xb is 



'^'^ From Eq. (12), there is 



Xb=Xa + 



f 

Jo 



dz 



(14) 
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where h is the medium thickness of FIG. 1 and FIG. 2. 
By substituting Eqs. (9) and (14)into (10), and using the 
equality 



n(0) sin6it^ = n{h) sinOl, 



we have 
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(15) 
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FIG. 3: The structure {BA)^ of the general function photoni( 
crystals. 



where 



6b = — nfe(6)(cos 6 + sin I 
c Jo 



dz 



I nliz) _ ^ 



= ), (17 



and similarly 



(18 



Substituting Eqs. (16) and (18) into (7) and (8), anu 

using H = ./^nE, we obtain 
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FIG. 4: The two periods transmission figure of light in general 
function photonic crystals. 



Ei 
Hi 



= Me 



En 
Hii 



(19) 



where 



Mr = 



cos^b 

-mfc(O)^/-^ cos 6l sin 6h 



i sin Sb 



Mb (0) cos cos (St IvJ^^i 

rib (^) cos 0^ 



The Eq. (20) is the transfer matrix M in the medium of 
FIG. 1 and FIG. 2. By refraction law, we can obtain 



smff = 



no 
n(0) 



sin 6^, cos Of 




n2(0) 



sin'^f, (21) 



where uq is air refractive index, and n(0) — n{z)\z=o- 
Using Eqs. (15) and (21), wo can calculate cos9(. 

4. The transmissivity and light field distribu- 
tion of one-dimensional general function photonic 
crystals 

In section 3, we obtain the M matrix of the half pe- 
riod. We know that the conventional photonic crystals 
is constituted by two different refractive index medium, 
and the refractive indexes are not continuous on the in- 
terface of the two mediums. We could devise the one- 
dimensional general function photonic crystals structure 
as follows: in the first half period, the refractive index 
distributing function of medium B is nb{z). and in the 
second half period, the refractive index distributing func- 
tion of medium A is na{z), corresponding thicknesses are 
b and a, respectively. Their refractive indexes satisfy con- 
dition ni,(b) 7^ na(0), their structure are shown in FIG. 
3, and FIG. 4. The Eq. (20) is the half period transfer 



matrix of medium B. Obviously, the half period transfer 
matrix of medium A is 



cos 6a 



Ma = 
where 



/lo „„(a) cose^J 

Ua (0) cos cos 5a 



-ina{Q)^ —COS Ol sin 6 a — ■ — , . ,,u 



U) 



6a = —na{a)[cosdi -a 
c 



■ sm 



'jii 



dz 



and 



cos el' 



sin 



"0 sin^ 0° 



3// 



riQ 
na{a) 



smt 



cos 



ill 



nl{a) 



■ sin^ e". 



In one period, the transfer matrix M is 
M = Mb ■ Ma 

cos 6b 



rib (b) cos Of 



-inb(0),f^ COS ehin 6b T^" 



-i sin Sa 



IJ-o 
COS^a 



-ina{Q).[^cosel'sin6a t c°^^. 

^ ' y A»o raa(a)cos«-^ 



^na(a) COS01' 



(23) 



(24) 



(25) 



(26) 



(27) 
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FIG. 5: The refractive index of the general functions in a 
period. 



FIG. 6: The relation between transmissivity and frequency 
corresponding to the general function photonic crystals. 



The defect layer's refractive index is constant n^, its 
transfer matrix is 



cos Sd — — sin Sd 



-irjd sin 5d cos 6d 



(28) 



where m = J j^n-d, 6. 



^Udd. 

The form of the GFPCs transfer matrix M is more com- 
plex than the conventional PCs. The angle 6(, 6f , Oj.^ 
and Oj^ are shown in Fig. 4. The characteristic equation 
of GFPCs is 



El 
Hi 



M1M2 ■■ ■ Md - ■ ■ Mn 



En+1 



MbMaMbMa ■■■Md--- MbMa 



M 



f E]si+i 
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N+l 
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.(29) 



Where N is the period number. With the transfer 

matrix M (Eq. (28)), wc can obtain the transmission 
and reflection coefficient t and t, and the transmissivity 
and reflectivity T and R, they are 



t = 



2% 



EtN+l _ 

En Arjo + Brjom+i + C + Dtjn+i ' 



r = 



Eri _ Atjo + BrjpTjN+i - C - Drjp 
En Ario + BrjoriN+i + C + Drjo ' 



and 



R 



t-t* 



r ■ r . 



(30) 



(31) 



(32) 



(33) 



In the following, we 



Where % = Vn+i = 

give the electric field distribution of light in the one- 
dimensional GFPCs. The propagation figure of light in 
one-dimensional GFPCs is shown in FIG. 9. From Eq. 
(28), we have 



El 
Hi 



■- Mi{di)M2{d2) ■ - ■ Md{d) ■ ■ ■ Miv-i(djv-i) 

M;v(A.;v) ( t J • • • t t a'^I ) (34) 
"'\HN(di+d2 h rfjv-i + A2;jv) / ^ ^ 



where di and ^2 are the thickness of first and second 
period, respectively, A.zn is the propagation distance of 
light in the N-th period, Ei and Hi are the intensity of 
incident electric field and magnetic field, and £'jv(rfi + 

d2 h dN~i + Azn) and HN{di + d2- ■ ■ + rfw-i + Azn) 

are the intensity of the N-th period electric field and mag- 
netic field. The Eq. (34) can be written as 

/ EN{di+d2 \-dN-i + AzN)\_,^-i,.^ X 

V H^idi + rf2- • • + dN-i + Azn) j " ^^^^^ 

M^l,{dN-i) ■ ■ ■ +M-\d) ■ ■ ■ M^\d2)M^\di) 

( A{Azn) B{Azn) \f Ei\ , . 

[ C{Azn) D{Azn) )\Hi)^ ^•^^> 

the electric field Ei and magnetic field Hi can be written 
as 



Ei = Eii+Eri = {l + r)Eii, 



Hi = Hiicosef - Hricosef 
^cos9^{l-r)Eii. 



(36) 



(37) 



From Eqs. (34)-(36), we can obtain the ratio of the elec- 
tric field EN{di+d2 \-dN-i + AzN) within the GFPCs 

to the incident electric field En, it is 

Eil 

= \A{Azn){1 + r) + B{AzN).f^cose^{l - r)|^8) 

V Ato 

5. Numerical result 

In this section, we report our numerical results of trans- 

missivity and light field distribution. We consider refrac- 
tive indexes of the linearity functions in a period, it is 

n(,(z) = n6(0) + ^^^fc^z, 0<z<b, (39) 



n4z) = na{0) + ^^^^^^^^^z, 0<z<a, (40) 

Eqs. (38) and (39) arc the line refractive indexes distri- 
bution functions of two half period mediums B and A. 
When the endpoint values n(,(0), n(,(6), na{0) and «„(«) 
are all given, the line refractive index hmctions nb{z) and 
na{z) arc ascertained. The main parameters are: the half 
period thickness b and o, the starting point refractive in- 
dexes 716(0) and na(0), and end point refractive indexes 
rib (6) and na{a), the optical thickness of the two medi- 
ums are equal, i.e., ni,{0)b = na(0)a, the incident angle 
9^ = 0, the center frequency ujo = 1.215 x 10^°Hz, the 
center wave length Aq = the thickness b = 280nm, 
a = 165nm and the periocl number N = 16. we take 
nb(0) = 1.38, ni){b) = 1.9 for the medium B, and 
'^a(O) = 2.35, na{a) = 2.6 for the medium A, which are 
the up line function of refractive indexes, it is shown in 
FIG. 5. By the refractive indexes function, we can cal- 
culate the transmissivity, we obtain the transmissivity 
distribution in FIG. 6. 

In FIG.7, we take nh(0) = nb{b) = 1.38, njO) = 
na{a) — 2.35, i.e., the transmissivity of conventional pho- 
tonic crystals. Compare FIG. 6 with FIG.7, it can be 
found the results: (1) when the line function of refrac- 
tive indexes is up, i.e., the GFPCs, the transmissivity can 
be far larger than 1 {T maximum is about 850), while 
the maximum of transmissivity in conventional photonic 
crystals is 1; (2) The number of band gaps in GFPCs are 
more than the conventional PCs. 

In FIGs.8-10, we main discuss the relation between 
transmissivity and wave length corresponding to the GF- 
PCs with defect layer, there is Udd = Base on the 
structure of GFPCs (BA)^^, we inset the defect layer at 
different position, which are shown in FIG.8(a-c). The 
structures are: {a){BAfD{BAf , {h){BAy°D{BAf , 
{c){BAy'^D{BA)'^, where = 2. We can see that when 
the position of defect layer is changed, the defect model 
become small. In FIG. 9, we compare the transmissivity 
with different structure as {BA)^D{BA)^. In FIG.9(a- 
c), the half period number are: N = 6, N = 7 and 
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FIG. 7: The relation between transmissivity and frequency 
corresponding to the conventional photonic crystals. 




FIG. 8: Comparing the transmissivity of the GFPCs with 
different position of defect layer. 



= 8, where Ud = 2. We can obtain the results: As the 
number of half period increase, for example, iV = 6, 
the maximum of defect model achieve 110 (FIG. 9(a)). 
When N = 7, the intensity of defect model is about 
250 (FIG.9(b)). When N increases up to 8 (FIG.9(c)), 
the maximum of defect model nearly 400, i.e., as the 
half period number N increase, the maximum of defect 
model increase, the defect model's width of half height 
become narrow and the transmissivity of the GFPCs also 
increase. 

FIG. 10 shows the transmissivity of the GFPCs 
{{BAfD{BAf) with different refractive indexes of de- 
fect layer. Prom FIG.lO(a-c), na are taken as: 1.5, 2 and 
5, respectively. Here, we notice that as the refractive 
indexes of defect layer increase, the intensity of defect 
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FIG. 9: Comparing the transmissivity of the GFPCs 
{{BA)^D{BA)'^) with different haff period numbers (a) N=6 
(b) N=7 (c) N=8. 



FIG. 11: The light distribution in the conventional PCs. (a) 

without the defect layer((S^)^^), (b) with the defect layer 
{{BAfD{BAf). The bold line is the field distribution of 
defect layer. 
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FIG. 10: Comparing the transmissivity of the GFPCs 
{{BA)^D{BA)^) with different refractive indexes of defect 
layer. 



FIG. 12: The light distribution in the GFPCs. The structures 
are: (a) (b) {BAfD{BAf, (c) {BAy''D{BA)\ The 

bold line is the field distribution of defect layer. 



model decrease, and the position of defect model become 
blue shift. 

FIGs. 11-14 are the distribution of electric field. The 
transverse axis z is propagation distance, and the longi- 
tudinal axis is the ratio of field intensity E and incidence 
field intensity En square, i.e., lE/Enl"^. 

FIG. 11 is the distribution of electric field in conven- 
tional PCs. The structure of FIG.ll(a-b) are (BA)"' 
and {BA)^D{BA)^. The main parameters arc the same 
as FIG. 7, and the = 2, n^d = It can be found that 
the electric field was enhanced obviously when inserted 
the defect layer. 

In the following, the main parameters are the same 
as FIG.6. In FIG.12, we study the effect of differ- 
ent structure of GFPCs on the distribution of light 
field. For FIG.12(a), the structure is (BA)^^, FIG. 12(b) 



and FIG.12(c), the structure arc {BAfD{BAf and 
{BAy^D{BA)^, where na = 2, Udd = ^. Comparing 
FIG.12(a) with FIG.12(b), we found that when the de- 
fect layer is located in the middle of GFPCs, the electric 
field was weaken. While in FIG. 12(c), the electric field 
was enhance. They are shown that both the defect layer 
and its' position have the effects on the distribution of 
light field, and we can find the defect layer made the elec- 
tric field local enhanced for the conventional PCs. While 
the defect layer made the electric field whole enhanced 
or decreased for the GFPCs. 

In FIG. 13, we shall consider the effect of half period 
number {{BA}^ D{BA)^) on the electric field. Taken 
AT = 4, AT = 6 and A/' = 8 corresponding to FIG.13(a- 
c) . The results shown that when N increase, the relative 
intensity of electric field heighten. 
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FIG. 14 show the effect of different refractive indexes 
of defect layer (n^;) on the electric field. The structure 
is {BA)^DiBAf. In FIG.14(a-c), Ud are equal to 1.5, 2 
and 3, respectively. As we can see in FIG.14(a-c), when 
rid increase, the relative intensity of electric field was en- 
hanced. 



6. Conclusion 
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FIG. 13: The light field distribution of the GFPCs 
{{B A)'^D{B A f)with different periodicity, (a) N=4, (b) 
N=6, (c) N=8. The bold line is the field distribution of defect 
layer. 




FIG. 14: The light field distribution of the GFPCs 

{(BA)^ D(BA)'^)with different refractive indexes of defect 
layer. The bold line is the field distribution of defect layer. 



In summary, We have theoretically investigated a new 
general function photonic crystals (GFPCs) with defect 
layer. Based on Ferinat principle;, we achieve the mo- 
tion equations of light in one-dimensional general func- 
tion photonic crystals, and calculate its transfer matrix. 
We choose the line refractive index function for two medi- 
ums A and B, and obtain some results: (1) When the po- 
sition of defect layer move behind in the GFPCs, or the 
refractive indexes of defect layer increase, the transmis- 
sion intensity maximum of defect model decreases. (2) 
When the period number of GFPCs increase, the trans- 
mission intensity of defect model increase, and the defect 
model's width of half height become narrow. (3) Both 
the defect layer and it's position have the effects on the 
electric field. (4) In the structure (BA)^ D{BA)^ , when 
the period number N increase, the relative intensity of 
electric field increase. (5) When the refractive indexes of 
defect layer (ud) increase, the relative intensity of electric 
field was enhanced. Since the GFPCs has new character 
different from the conventional PCs, it should be helpful 
in the design of photonic crystals. 
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